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Reorganization asymmetry of electron transfer in ferroelectric media 
and principles of artificial photosynthesis 

Dmitry V. Matyusho\Q 

Department of Chemistry and Biochemistry and the Center for the Early Events in Photosynthesis, 
Arizona State University, P0 Box 871604, Tempe, AZ 85287-1604 
(Dated: February 6, 2008) 

This study considers electronic transitions within donor-acceptor complexes dissolved in media 
with macroscopic polarization. The change of the polarizability of the donor-acceptor complex 
in the course of electronic transition couples to the reaction field of the polar environment and the 
electric field created by the macroscopic polarization. An analytical theory developed to describe this 
situation predicts a significant asymmetry of the reorganization energy between charge separation 
and charge recombination transitions. This result is proved by Monte Carlo simulations of a model 
polarizable diatomic dissolved in a ferroelectric fluid of soft dipolar spheres. The ratio of the 
reorganization energies for the forward and backward reactions up to a factor of 25 is obtained 
in the simulations. This result, as well as the effect of the macroscopic electric field, is discussed in 
application to the design of efficient photosynthetic devices. 



I. INTRODUCTION 

Electron transfer (ET) reaction is a basic elementary 
step in photoinduced charge separation occurring in both 
natural and artificial photosynthesis. The fundamental 
mechanism behind the conversion of light energy into the 
energy of a charge-separated state is illustrated in Figure 
n Optical excitation hv of the donor unit, D— >D*, within 
the donor-acceptor complex, D-A, lifts the energy level 
of the donor to create conditions for the photoinduced 
charge separation step, D*-A — > D + -A~. Charge separa- 
tion is generally an activated transition. The activation 
barrier, according to the Marcus-Hush theory?! can be 
obtained from the equilibrium free energy gap between 
the final and initial ET states, AFo, and the free energy 
A required to reorganize the nuclear degrees of freedom 
(reorganization free energy). When the equilibrium free 
energy of the acceptor is below the free energy of the 
donor by the amount of free energy A, the charge sepa- 
ration transition is activationless. Such an activationlcss 
transition is often fast and efficient, but it requires loosing 
the energy A from the photon energy hv. Increasing the 
energetic efficiency of photosynthesis therefore demands 
A to be as low as possible. It is often assumed that an 
important role played by the hydrophobic protein matrix 
in facilitating ET in natural photosynthesis is to screen 
highly polar aqueous environment and to reduce A. 

Once the charge-separated state has been created, it 
needs to be sufficiently long-lived to be used for energy 
storage. A mechanism often suggested to be of primary 
importance in natural systems is the use of a sequence of 
activationless ET steps to move the electron away from 
the primary donor. Each activationless step i lowers the 
system energy by A^ thus resulting in the overall low ener- 
getic efficiency of the photosynthetic unit. The reaction 
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competing with photoinduced charge separation is the 
normally highly exothermic return electron transfer to 
the ground D-A state. The classical reaction path then 
goes through the relatively high activation barrier in the 
inverted ET region predicted by the Marcus-Hush theory 
(point "C" in Figure^. 1 However, the system avoids this 
activated path by nearly activationless transfer to a vi- 
brationally excited state, D-A(v), of the electronically 
ground donor-acceptor complex (dashed line in Figure 
[TJ. The vibrational energy is subsequently released to 
heat through vibrational relaxation (v — > in Figure 

The design an optimization of a sequence of activation- 
less ET reactions requires a very precise molecular tun- 
ing, which is hard to achieve in synthetic systems despite 
some significant progress achieved in this field in recent 
yearsi^ It is therefore desirable to search for mechanisms 
of reducing the rate of return ET in molecular photosyn- 
thesis with the goal of achieving higher quantum yield 
for the charge-separated state. According to the current 
understanding of radiationless transitions in molecules, 4 
an efficient way to reduce the return rate would be to 
move the reaction into the normal region of ET. 

The electronic states D*-A and D-A are distinct and, 
in principle, the charge-recombination transition D + 
A~ — * D-A can be characterized by a pair of parabolas 
with the curvatures different from those of the charge- 
separation transition D*-A — > D + -A~. In the reactions 
classification adopted here we label transition D*-A — > 
D + -A~ as charge separation and transition D + -A~ — > 
D-A as charge recombination. Backward transitions at 
each step are not considered as separate steps in the re- 
action mechanism. When the photoexcitation energy is 
kept constant, shifting the recombination reaction into 
the normal region would require increasing the curva- 
ture of the charge-recombination parabolas (dash-dotted 
lines in Figure^), i.e. lowering the reorganization energy. 
The inverted-region activated state C then shifts to the 
normal-region activated state C in Figure ^ 

It is easy to realize the pitfall of the picture shown 
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in Figure ^ The separation of the minima of two dash- 
dotted parabolas must be equal, in the Marcus-Hush the- 
ory, to twice the reorganization energy, which is clearly 
violated once the curvature is increased without corre- 
sponding shift of the parabolas. Therefore, the goal of 
bringing the recombination reaction to the normal region 
can be realized only within models extending beyond the 
Marcus-Hush picture of equal-curvature parabolas. On 
needs flexibility, built into the model, that would allow 
decoupling of the Stokes shift from the curvatures. The 
use of parabolic free energy surfaces with different curva- 
tures, as was suggested by Kakitani and Mataga, 5 is pro- 
hibited by the requirement of energy conservation. When 
the energy gap between the donor and acceptor energy 
levels is taken for the reaction coordinate, the free en- 
ergy surfaces of the initial, Fi(X), and final, F<z{X), ET 
states are connected by the linear relation established by 
Warshel 6 - and Tachiyai 

F 2 (X) = F^X) + X (1) 

The parabolic surfaces with different curvatures clearly 
violate this requirement and, therefore, cannot be used 
for the modeling of ET reactions. The problem can be re- 
solved within a three-parameter model of ET free energy 
surfaces 8 which allows different reorganization energies 
and, at the same time, does not violate eq^ The free 
energy surfaces Fi(X) are then necessarily non-parabolic. 

Once different reorganization energies are allowed 
within a theory of free energy surfaces, one needs to ad- 
dress the following major question: What would be a 
mechanism resulting in dramatic changes of reorganiza- 
tion energies between charge separation and charge re- 
combination reactions? A question relevant to studies 
of natural photosynthesis is the following: Is the role of 
the protein matrix can be reduced to providing a low 
dielectric constant or there are some other properties of 
proteins beneficial for high efficiency of natural photosyn- 
thesis? Note that shifting the transition point from C to 
C requires a very significant change in the reorganization 
energy (four times decrease in Figure ^ . The possibil- 
ity of different equilibrium-point curvatures of the two 
ET free energy surfaces has been actively pursued in the 
last decades by computer simulations searching for the 
effects of non-linear solvation on ETi 6 ! 9 ; 10 ! 11 ! 12 ! 13 How- 
ever, essentially in all reports, the difference in reorgani- 
zation energies between the initial and final ET states has 
never come close to the magnitude that would dramati- 
cally change the mechanism of charge recombination. 

In this paper, we address the problem of reorganiza- 
tion asymmetry by considering ET reactions in polar- 
izable donor-acceptor complexes immersed in polariza- 
tion anisotropic media. A part of the motivation for 
this approach comes from studies of bacterial photosyn- 
thesis. The protein matrix in natural bacterial reac- 
tion centers is highly anisotropic with the electric field 
amounting ~ 10 7 V/cmi^ On the other hand, the pri- 
mary pair is highly polarizable with the polarizability 
change upon photoexcitation about 800-1100 A 3 ii^ Since 
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FIG. 1: Energetics of photoinduced electron transfer involved 
in natural and artificial photosynthesis. The charge-separated 
state D + -A~ is placed below the photo-excited state D*-A 
by the reorganization energy A to ensure activationless tran- 
sition. The dashed line indicates the vibrationally excited 
electronic ground state for which the recombination transi- 
tion D + -A~ — >D-A does not require an activation barrier. 
The dash-dotted lines mark the free energy surfaces of states 
D + -A~ and D-A with the reorganization energy four times 
lower than the value of A used to draw the free energy sur- 
faces shown by the solid lines. C and C' indicate the classical 
transition states. 



bacteriochlorophyll cofactors are much less polarizable, 16 
primary charge separation occurs with a large negative 
polarizability change. On a more fundamental level, po- 
larizability change in the course of ET leads to very sig- 
nificant reorganization asymmetry, as follows from both 
analytical theories^ 7 - and computer experiment The 
asymmetry arising from the coupling of changing polar- 
izability to the solvent recation field is by far the largest 
achieved so far in computer simulations thus promising 
the desired alteration of the energetics of photosynthetic 
reactions. 

In Section [HJ we present an analytical model of ET 
free energy surfaces when both the polarizability and the 
dipole moment of the donor-acceptor complex change in 
the course of ET. The donor-acceptor complex is im- 
mersed in a polar medium which, in addition to usual 
polar response, is characterized by a macroscopic elec- 
tric field. In order to mimic high electric fields present 
in molecular systems (thin films, proteins, etc.), we use a 
polar model fluid capable of producing ferroelectric order 
( Section ITTT|l . The ordered phase is obtained by Monte 
Carlo (MC) simulations. The possibility of ferroelectric 
order in bulk polar fluids is still actively debated in the 
Iiteraturej22i2ii2£ and it has been suggested that bound- 
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ary conditions employed in the simulation protocol sig- 
nificantly affect the possibility of creation of the macro- 
scopically polar phased For our current purpose, this 
subject is not relevant since we are mostly interested in 
reactions in condensed phases with macroscopic electric 
fields strong enough to influence properties of electronic 
transitions. Whether the macroscopic polar order is sta- 
bilized by surface effects or some other reasons is not 
central to our study. On the other hand, ferroelectrics 
and ordered polar films may be used as solvents for ET 
reactions. The present analysis then directly applies to 
such systems. We return to discussing the role of reor- 
ganization assymetry in improving the efficiency of ph- 
tosynthesis in Section HVI 



where a„ is the linear response function of the solvent 
nuclear polarization. In the following, for simplicity, we 
will assume collinear dipole moments and the polar- 
izability changing its value only along the direction of 
the dipole moment. We will also assume that the so- 
lute polarizability increases for the 1 — > 2 transition, i.e., 
Act = ot2 — oti > 0. These approximations allow us to 
consider a„ as a scalar with the only non-zero projection 
along the solute dipole. We will also define the scalar F m 
as the projection of the external field F on the direction 
of the solute dipole moment. The consideration of a more 
general case does not present fundamental difficulties^ 



B. Free energy surfaces 



II. ENERGETICS OF ELECTRON TRANSFER 
A. Model 

Consider a donor-acceptor complex immersed in a 
medium with a nonzero macroscopic electric field F. The 
complex has the dipole moments moi and polarizabili- 
ties a.Qi in two electronic states, i = 1,2. These states 
are coupled to nuclear fluctuations in the solvent and are 
"dressed" with the field of the electronic solvent polariza- 
tion following adiabatically the changes in the charge dis- 
tribution within the complex. Because of the electronic 
polarization, the dipole moments and polarizabilities are 
different from their gas phase values moi and aoi*^ 

-l 



mi = [1 - 2a ai ■ a e • m 0i 



(2) 



and 



ol 1 = [1 - 2a 0l ■ a e ] • a Qi 



(3) 

Here, a e is the linear response function such that the 
free energy of solvation of dipole mo by the electronic 
solvent polarization (subscript "e") is mo • a e ■ mo. Once 
the electronic degrees of freedom have been adiabatically 
eliminated, the electronic energy levels of the charge- 
transfer complex are affected by the microscopic nuclear 
field R n (nuclear reaction field, subscript "n") and the 
macroscopic field F of the solvent *i£ 

Ek [Rn] = If - ' (R« + P) ~ \ (R« + F) ■ a* ■ (R„ + F) 

(4) 

Here are the energies of the ET complex which in- 
clude the gas-phase energies and free energies of solvation 
by the solvent electronic polarization expressed as a sum 
of induction and dispersion solvation components. 
The system Hamiltonian 

flj[Rri] = Ei[R, n ] + i/e[R„] (5) 

is the sum of energies Ei [R„] and the solvent bath Hamil- 
tonian Hb- The Gaussian (linear response) approxima- 
tion is adopted for the latter 
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~Rr, 



Rk 



(G) 



The classical Hamiltonian Hi[TL n ] in eqJSJcan be used 
to build the free energy surfaces of ET along the reac- 
tion coordinate associated with the fluctuating donor- 
acceptor energy gap required to reorganize the nuclear 
degrees of freedom 



x = AH[n n ] = H 2 [n n ] - flifR*] 



(7) 



The free energy surfaces for the initial (i = 1) and final 
(i = 2) ET states are obtained by constrained integration 
over the nuclear reaction field 



e -m(x) = A 



J 6(X-AH[K n })e^ H ^dR n , 



(8) 



where A is used to account for the units of the field R„ 
and /3 = l/k B T. 




FIG. 2: Free energy surfaces of ET Fi(X) for the initial (i — 1, 
"1") and final (i = 2, "2") ET states. Shown in the plot are 
the free energy minima Xoi, the free energy gap AFo = F02 — 
Fqi, and the upper boundary Xo for the reaction coordinate 
X. 

The Gaussian integral in eq |H1 can be taken exactly 
with the result 



Fi(X) = -^X-p-Hn 



sinhxPO 



where 



X(X) 



(9) 
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The parameters in eas 1^1 and 1101 are related to the prop- 
erties of the polarizable donor-acceptor complex and the 
polar solvent by the following set of equations. The pa- 
rameters Ki are 



Ki = (2a n f ni Aa)~ 



(11) 



The parameter f n $ in ea llll is responsible for the enhance- 
ment of the effective solute dipole by the interaction of 
the solute polarizability with the nuclear reaction field 
(cf. to cqs |21 and 



fni = (1 - 2a„a 4 ) 1 



(12) 



Further. 



Xq = - 7- " ^F m - lAaFl + ^ ^Aa ^ 

(13) 

is the boundary of the fluctuation band of the reaction 
coordinate X (Figure [3J|. The requirement X < Xq im- 
plicit in the definition of xPO is eq E3 docs not allow 
the energy gap fluctuations to exceed Xq. Finally, Ai in 
eg HUl is the solvent reorganization energy 

A l = a n f ni [Am + Aa(R t + F)] 2 (14) 

where Am = m.2 — mi and reaction field Ri is 

Ri = 2a n f ni mi (15) 

The function x(^0 m eqEldoes not carry index speci- 
fying the ET state because of the relation connecting the 
parameters of the two ET states: 

KfXi = nl\ 2 (16) 

In addition, the parameters K\ and k 2 are related by the 
additivity constant 



K\ - K2 = 1 



(17) 



The free energy surfaces defined in the range X < Xq 
are non-parabolic, as is depicted in Figures|21and|3J Each 
of them passes through the minimum at X — X$i defined 
by the condition 



dFi(X)/dX 







(18) 



X=X 0i 



This equation can be reduced to the algebraic equation 
for x{X) which always has a non-zero solution 



1,,2 



X coth X -l = (2/3<A 4 )" i X 



(19) 



The solution is simplified when \oi = x(^oi) 1. One 
then gets 



Xoi — X — KiXi 
and the free energy at the minimum is 



Fen — —KiXo — \i 



(20) 
(21) 



In this limit, the second derivatives taken at the posi- 
tion of the minima define the two solvent reorganization 
energies, as is the case in standard formulations of ET 
theories 



d 2 F(X) 



dX 2 



X=Xn 



1 

2A~ 



(22) 



Note that from eq |^ the fluctuation boundary Xq is 
related to the free energy gap AFo = F02 — F01 and two 
reorganization energies by the relation 



Xq = AF + - 



(23) 



Taken together, eqs EH23 provide an exact model for 
the free energy surfaces of ET (Figure^ based on three 
thermodynamic parameters, Ai, A2, and AFq. The 
present model thus extends the two-parameter Marcus- 
Hush theory, based on the assumption Ai = A2J* to three- 
parameters space. In view of the connection between m 
and Ki (eq 1171) , one of them can be considered as the 
non-parabolicity parameter: 



K = K\ = \\ — \ — 




(24) 



In terms of the solute polarizability and the solvent nu- 
clear response function, the non-parabolicity parameter 
becomes 



1 - 2a n a 1 
2a„Aa 



(25) 



Because of the choice Act > adopted here, Ai < A2 and 

Ki > 0. 

It is easy to prove that Fi(X) from eq[!|]obey the en- 
ergy conservation requirement given by eq^ Note that 
the present model exhausts the space of thermodynamic 
parameters available for the description of ET reactions. 
Only three thermodynamic parameters, free energy gap 
between the minima and second cumulants of the fluctu- 
ations around the minima, are allowed by the two-state 
nature of the problem. Any extension of the theory for 
ET free energy surfaces beyond the present level will re- 
quire non-equilibrium parameters to be involved. 

The requirement of thermodynamic stability of the po- 
larization fluctuations limits the range of possible values 
of the reaction coordinate, X < Xq. The reaction coor- 
dinate boundary Xq is shown by the dash-dotted line in 
Figure |3 There is no singularity of Fi(X) at X = Xq 
and the free energies Fi(X) are infinite at X > Xq. For 
reaction coordinates sufficiently far from Xq, the free en- 
ergy surfaces can be conveniently re-written in the form 



F(X) - F 0i 



X 0i + Ki Xi - X - \J K t \ t (26) 



When the reorganization energies A^ are close to each 
other, the non-parabolicity parameter k in eq 1241 tends 
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FIG. 3: Free energy surfaces F;(A") at Ai = 1 eV and A2 = 2 
eV. The free energy gap AFo is: —2 eV and —7 eV in (a) and 
1 eV and 4 eV in (b). The free energy surfaces corresponding 
to AFo — — 7 eV in (a) do not have the classical crossing point 
X = indicated by the dashed lines in both panels. 



to infinity and one obtains the standard Marcus-Hush 
parabolas by expanding ea 1261 in 1/k. From ea 1261 the 
activation energy of ET is 



AEi = Fi(0) - Fm 



\J X 0i + KiXi — \J K. t Xi (27) 



C. Qualitative results 

The three-parameter model predicts some novel results 
regarding the dependence of the reaction rates on the 
free energy gap (energy gap law). In order to illustrate 
them, we will consider the transition from the activated 
normal region to the activated inverted region through 
the activationless transition for the forward reaction 1 — > 
2 and the backward reaction 2 — > 1 while maintaining 
Ai < A2 and the definition of the free energy gap as 
AF = F 02 — Fq\ (Figure El . Our analysis here assumes 
Xoi 1, which holds for most cases of interest. For the 
forward reaction 1 — * 2, the transition from the normal to 
the inverted region is marked by the equation Xq\ = 0, 
which corresponds to a line in the space of parameters 
A2/A1 > 1 and AF /\i (Figure E}: 



AFq/Ax = k(k + !)-(«- 1) 2 (A 2 /Ai) 



(28) 



The line shrinks into the point AFq/Xi = ±1 in 
the Marcus-Hush limit of equal reorganization energies 
(marked as "MH" in Figure 0J). 

Lowering the free energy gap in the inverted region for 
the forward reaction 1— > 2 leads to a new region of ET 



absent in the Marcus-Hush theory. When the boundary 
of the band of allowed energy gaps Xq crosses the transi- 
tion state X = 0, the classical crossing point of two free 
energy surfaces falls outside the range of allowed reaction 
coordinates. No crossing of free energy surfaces Fi(X) is 
then possible (AFo = — 7 eV in Figure Ola)), and the 
classical reaction channel is closed, AFi = 00. The re- 
action occurs only through vibrational excitations of the 
final ET state 2 effectively lowering the free energy gap. 
The transition to this new region, which may be called 
"quantum tunneling region" (AF X = 00 in Figure^Ja)), 
is marked by the line X = is the space of parameters 
A2/A1 and AFq/Ai. 
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FIG. 4: Normal and inverted region in the space of parameters 
AF0/A1 and A2/A1. Shown are the results for the forward 
reaction 1 — > 2 (a) and for the backward reaction 2 — > 1 
(b). The quantum tunneling region with AFi = 00 in (a) 
is separated from the inverted region by the condition Xo = 
0, which puts the classical crossing point outside the range 
of reaction coordinates X < Xo allowed by the condition of 
thermodynamic stability. The dots labeled as "MH" refer 
to the Marcus-Hush limit in which the transition from the 
normal to inverted region is given by the conditions: AFo / A = 
-1 for 1 -> 2 and AF /A = 1 for 2 -> 1. 

Due to the asymmetry of the free energy surfaces for 
A2/A1 > 1, the normal region spans different ranges of 
AF for positive and negative free energy gaps. Figure^] 
shows a broader normal-range ET for AFo > 0. This ob- 
servation suggests that the goal of bringing the exother- 
mic recombination reaction to the normal region (Figure 
is easier to achieve when the reorganization energy of 
the charge-separated state D + -A~ is much higher than 
the reorganization energy of the ground state D-A (see 
Discussion) . 

The overall classical energy gap law (no vibrational ex- 
citations) for the present model is illustrated in Figure 
With increasing the ratio of the reorganization energies 
A2/A1, the bell-shaped dependence of — AFi on AF be- 
comes increasingly shallow in its right wing with positive 
energy gaps AFo > 0, approaching the linear energy gap 
law AFi cx AFq. On the other hand, the approach of the 
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FIG. 5: Energy gap law for the reaction 1^2 (Ai < A2) at 
Ai = 0.2 eV and A 2 = 0.25 eV (dash-dotted line), A 2 = 0.4 
eV (dashed line), and A2 = 2 eV (solid line). The dotted 
lines separates the normal (N), inverted (I), and quantum 
tunneling (QT) regions for A2 = 2 eV. 



The Arrhenius factor for the rate constant can be ob- 
tained (within a factor) by putting X = in eq|2U Since 
the thermodynamic stability of the nuclear fluctuations 
requires X < Xo n , only terms corresponding to Xon > 
in eq|3n|will contribute to the sum in eq[2Hl This fact im- 
plies a modification of the standard picture of the FCWD 
made of a Poisson- weighted sum of vibronic transitions. 
Vibronic transitions with n < —Xo/(huj v ) will not con- 
tribute to the FCWD as is illustrated by n = and n = 3 
in Figure [SJ Only starting from vibrational excitation 
number making Xo n positive (n — 6 in Figure will a 
given vibronic transition contribute to the FCWD. This 
fact will result in a greater asymmetry of emission lines 
making the red-side wing of an optical band more shallow 
than the blue-side wing. 



fluctuation boundary Xq to the transition point X = 
squeezes the left wing narrowing the inverted region of 
ET. 
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FIG. 6: Inverted-region reaction 1 — > 2 with Ai = 0.2 eV, 
A2 = 1.0 eV, and AFq = —2 eV (shown by the vertical ar- 
row) . Classical transitions between vibrationally ground-state 
surfaces are forbidden by the condition Xq < 0. The dashed 
surfaces obtained for n = 3 do not contribute to the Franck- 
Condon weighted density of states since Xq + nhuj v < for 
n — 3. For n > 6, Xo + nfiuj v > (hcu v = 0.2 eV) and these 
vibronic transitions contribute to the overall density of states. 
The dash-dotted lines indicate the fluctuation boundary Xo n 
at n = and n = 3; X 0n = at n = 6. 

The quantum vibronic excitations of the donor- 
acceptor complex can be included in the standard way 
by considering the Franck-Condon weighted density of 
states (FCWD) as a Poisson-weighted sum over the vi- 
bronic excitations leading to n vibrational quanta in the 
final ET stated For instance, for the forward reaction one 
gets 

00 Oil 

FCWD(X;1^2) = e - s V * e-^i-W (29) 

Here, S = Xi/hu> v is the Huang-Rhys factor, u> v is the 
characteristic vibrational frequency, and F\ n (X) is ob- 
tained from eqlSJby replacing there Xq with 

Xo n = Xq + nhu> v (30) 



III. ELECTRON TRANSFER IN 
FERROELECTRICS 

Equation ^] for the solvent reorganization energies of 
a polarizable donor-acceptor complex is the central re- 
sult of the formal theory relevant to our discussion of 
reorganization anisotropy in ferroelectrics. Equation 1141 
predicts that the polarizability change Aa couples to 
the macroscopic and local (reaction) fields to modify the 
change in the solute dipole moment in the course of elec- 
tronic transition. Depending on the relative signs of Am 
and Aa(Ri + F m ), the solute polarizability may increase 
or decrease the solvent reorganization energy compared 
to the non-polarizable limit. The main question in this 
regard is whether the term Aa(Ri + F m ) can achieve a 
magnitude comparable to Am. Below, we address this 
question by MC simulations of a model donor-acceptor 
diatomic in a dipolar ferroelectric solvent. Here, we first 
provide some relevant estimates. 

The mean-field Weiss theory of spontaneous dipolar 
polarization 21 relates the macroscopic field in a sample 
of aligned dipoles m to their number density p as 

F = (4tt/3)pto (31) 

The term AaF m /m at p* = pa 3 = 0.7 is then ~ 3Aa/a 3 
when the solute dipole is aligned with the macroscopic 
field, a is the solvent diameter. The polarizability change 
may vary substantially between different donor-acceptor 
units, but the ground-state polarizability is close to Cq/16 
for many molecular systems (00 is the effective solute di- 
ameter). If Aa is of the same order of magnitude as the 
ground-state polarizability, 3Aa/cr 3 scales as (a /a) 3 . 
For example, for primary charge separation in Rhodobac- 
ter sphaeroids, Am ~ 53 Di^i If we accept a — 2.87 A 
and m = 1.83 D for water and Aa ~ 800 Af£ then 
Am/m ~ 29 and AaF m /m ~ 100 thus resulting in 
comparable order-of-magnitude contributions from the 
change in the permanent charges and from the interaction 
of the polarizability change with the non-homogeneous 
electric field. In addition, from the estimate of the lo- 
cal electric field of the protein matrix at the position of 
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the primary pair— F ~ 10 7 V/cm, the induced dipole 
AaF ~ 27 D is also comparable to Am ~ 53 D. 

A. Simulations 

A fluid of soft spheres is known to spontaneously form a 
ferroelectric liquid phase when conducting boundary con- 
ditions are employed in the simulation protocoL 20 This 
system was used here to model a disordered phase with a 
molecular-scale macroscopic polarization coupled to the 
electronic states of the donor-acceptor complex. Two 
sets of simulations have been carried out. The fluid of 
soft dipolar spheres was simulated in the first set in or- 
der to establish the range of parameters for which fer- 
roelectric phase can be detected. This was followed by 
the second set of simulations in which polarizable donor- 
acceptor complex was dissolved in the ferroelectric liquid. 

MC simulations of the pure solvent employed NVT en- 
semble of N = 500 particles in a cubic simulation cell at 
the reduced density of p* = pa 3 = 0.7. The molecules 
are interacting with the potential 

w(12) =4e((i/ri 2 ) 12 -mi -Ti2 -ma (32) 

where e is the repulsion energy, a is the diameter, nij 
is the dipole moment with the magnitude to, T12 = 
ViVar^ 1 is the dipolar tensor, and ri 2 = |i"i — r 2 |. Sim- 
ulations were done at fie — 1.35 and varying reduced 
dipole (to*) 2 = j3m 2 /a 3 . Periodic boundary conditions 
and reaction-field correction for dipolar interactions with 
conducting boundary conditions 2 ^ have been employed. 
The length of simulations was 10 6 cycles long far from 
the transition point and up to 10 7 cycles long close to 
the transition to ferroelectric phase. 

Transition to ferroelectric phase was monitored by cal- 
culating the first-order, S\ , and second-order, S2 , param- 
eters (FigureUfa)). The first-order (ferroelectric) param- 
eter quantifies the spontaneous polarization 

S 1 ! = M/Nm (33) 

where M is the total dipole moment of the solvent. 
The second-order (nematic) parameter is defined as the 
largest eigenvalue of the ordering matri?* 2 ^ 

Q = (2A)- 1 ^(3e J e J -l) (34) 
3 

where &j = xa.j/m. 

Both order parameters change smoothly with to* and 
do not allow a reliable identification of the transition 
point— Susceptibilities, which are expected to show 
sharp spikes at the points of phase transition^ are better 
indicators. Indeed, the dielectric susceptibility 

Xp = ((3/V)((SM) 2 ) (35) 

shows the first peak at (m* F ) 2 — 6.9 and the second peak 
at (m* c ) 2 = 7.5 (Figure □»), where V in eqESJis the 
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FIG. 7: First (Si, circles) and second (S2, squares) order 
parameters of the fluid of dipolar soft spheres vs (3m 2 /a 3 (a). 
Also shown are the dipolar susceptibility \p (b) and constant- 
volume heat capacity cv (c). The dotted lines connect the 
simulation points. The dash-dotted lines in (b) indicate the 
points of phase transition to ferroelectric fluid (nip) and fee 
crystal (m*?). 

solvent volume. The heat capacity 

c v =S/2 + (3 2 ((SE) 2 )/N (36) 

on the contrary, shows only a little bump at m* F and 
a strong peak at m* c (Figure 0c)), where E in eci I5B1 
is the total energy of the fluid. The order parameters 
also show weak discontinuities at vn* c . The examination 
of the density structure factors reveals that the system 
is in the liquid state below m* c and crystallizes into an 
fee lattice at this value of the reduced dipole. The point 
m F corresponds to the transition to a ferroelectric liquid. 
Note that the value of m* F obtained here is very close 
to that predicted by a linear extrapolation of the phase 
transition line recently reported by Weis 2 ^ for dipolar 
hard spheres. 

The ferroelectric liquid at (to*) 2 = 7.0, p* — 0.7, and 
(3e = 1.35 was used to study the dependence of the sol- 
vent reorganization energy on solute polarizability. In 
these simulations, a donor-acceptor diatomic made of 
two hard fused spheres of diameters uo/a — 1.5 with 
the center-to-center separation d/a — 0.6 was inserted in 
the center of a cubic simulation cell containing N = 500 
soft dipolar particles at p* = 0.7. Two opposite charges 
<Zo = ftllol/v = 10 where placed at the centers of two 
spheres in the charge-separated state thus producing the 
dipole moment TO02 = d\qo\- The induced point dipole 
moment po; = eoaoi-Ri, caused by the dipolar polariz- 
ability aoi and reaction field Ri, is placed at the mid- 
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FIG. 8: Charges and polarizabilities of the donor-acceptor 
complex in three states involved in photoinduced charge sep- 
aration (CS) and charge recombination (CR). 



point of the line connecting the two centers of the di- 
atomic, poi is aligned along the unit vector eo along the 
same line. The solute induced dipole was equilibrated to 
each instantaneous configuration of the solvent by itera- 
tion algorithm (more details on the simulation protocol 
are given in ref ll9T) . Ewald sums with conducting bound- 
ary conditions^ were used to calculate the solute-solvent 
electrostatic interactions. 

Some complications arise from the fact that the di- 
rector, d = M/7Vm, fluctuates in the simulations. In 
real liquid crystals, fluctuations of the macroscopic direc- 
tor occur on the time-scale of microseconds to seconds^ 
much slower than orientational motions of molecular so- 
lutes. The director, therefore, does not fluctuate on the 
time-scale of the reaction. In order to account for this hi- 
erarchy of relaxation times, the simulation protocol was 
set up to keep the orientation of the donor-acceptor com- 
plex fixed relative to the director. The initial configura- 
tion was created by inserting a sphere of diameter Co 
into the soft dipolar solvent followed by a short, 5 x 10 3 
cycles, equilibration run designed to establish the direc- 
tor d. The two spheres making the diatomic were then 
pulled apart at a given angle to the director and this 
angle was adjusted after each cycle over the N solvent 
molecules. The data were collected from simulation runs 
of the length 6x 10 5 — 1.2x 10 6 cycles. In the present sim- 
ulations, the long axes of the diatomic was always aligned 
with the director in such a way that the charge-transfer 
dipole Amo is always parallel to the macroscopic field F. 



B. Energetics of photoinduced electron transfer 

We have followed the basic design of photoinduced ET 
outlined in Figure ^ and detailed in terms of charges and 
polarizabilities in Figure[S] It is assumed that the donor- 




0.01 



FIG. 9: The linear response parameter £ from eg I4UI (a) and 
the distribution function of the ferroelectric order parameter 
(b) obtained from MC simulations at (m*) 2 = 7.0, f3e = 1.35, 
and p* = 0.7. The dotted line in (a) connects the simulation 
points. 



acceptor complex has no appreciable dipole moment and 
polarizability in its ground state (qo = 0, ot g = 0). Pho- 
toexcitation of the donor to state D*-A lifts the polar- 
izability to aoi > but does not substantially change 
the charge distribution, qo = 0. Both the charge distri- 
bution and polarizability change upon charge separation 
resulting in D + -A~ state (qo > 0, a 02 > 0). Because of 
the coupling of the polarizability to the reaction field Ri 
and the macroscopic electric field F m (ea 1140. the reor- 
ganization energies will differ for forward and backward 
electronic transitions for both the charge separation and 
charge recombination steps of the reaction mechanism. 
Therefore, we need four reorganization energies, Ai and 
A2 for charge separation and A3 and A4 for charge recom- 
bination (Figure EJ. 



TABLE I: Energetic parameters (eV) of the charge separation 
reaction (Figure |SJ in which aoi is varied and aro/o" 3 = 0.2 



is kept constant; (3qo/a 
QQl/cr 3 



A 



10 in the charge-separated state. 



0.2 
0.3 
0.4 
0.5 
0.6 
0.7 
0.8 



0.324 0.456 

0.344 0.371 

0.373 0.295 

0.477 0.229 

0.587 0.171 

0.761 0.123 

1.314 0.083 



0.817 
0.750 
0.706 
0.704 
0.665 
0.612 
0.645 



0.770 
0.715 
0.665 
0.666 
0.647 
0.640 
0.733 



"Calculated from the simulation data assuming (3 = 40 eV 1 . 
''From th e sim ulation data. 
c From eqlSSl 



The reorganization energies were calculated from MC 
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configurations as second cumulants of the solute-solvent 
interaction potential^ The polarizability of the photoex- 
cited state cxqi was varied at constant ace/c 3 = 0.2 for 
charge separation and the polarizability of the charge- 
transfer state ao2 was varied at a g — for charge re- 
combination. Tables [I] and [H] and Figure EH summa- 
rize the forward and backward reorganization energies 
for charge-separation and charge-recombination reactions 
along with the Stokes shift 



AXq — Xqi — Xq2 

From ea 1201 the Stokes shift becomes 

AX = k(X 2 - Ai) - A 2 



(37) 



(38) 



The Stokes shift in eqEZlis equivalent to the difference of 
absorption and emission maxima due to eq^ Equation 
1381 also applies to the Stokes shift defined in terms of two 
first spectral moments 



AX = (X) 1 -(X) 2 
in the limit myPXi ^> 1. 



(39) 
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FIG. 10: Reorganization energies and the Stokes shifts for 
the charge separation (CS) and charge recombination (CR) 
reactions (Figure[SJ vs polarizability of the photoexcited state 
aoi (CS) and the polarizability of the charge-separated state 
a 02 (CR); 9o * = 10, a g = 0. 

The direct comparison of the results of simulations to 
the three-parameter model discussed in Sectionals com- 
plicated by the narrow range of solvent dipoles for which 
ferroelectric phase can be detected. The close proximity 



of two phase transition points makes the statistics of po- 
larization fluctuations non-Gaussian, in contrast to the 
Gaussian approximation adopted in eqEJ There are sev- 
eral indications of significant deviations from the Gaus- 
sian statistics. The ratio 

Z = -(v ss )/f3((5v ss ) 2 ) (40) 

is equal to one (£ = 1) when polarization fluctuations are 
Gaussian (linear response approximation)^ Here, (v ss ) 
is the average electrostatic interaction energy of a liq- 
uid particle with the rest N — 1 particles in the liquid 
and ((Sv ss ) 2 } is the variance of this electrostatic interac- 
tion. As is seen from Figure^a), this parameter is higher 
than one in the paraelectric phase and drops sharply at 
the transition to the ferroelectric phase. In addition, the 
distribution of macroscopic polarization M in the simu- 
lation box is non-Gaussian. A shoulder seen in Figure 
HJb) points to the importance of a oc M 4 term in the po- 
larization functional such as that present in the Landau 
theory of phase transitions i2£ 

In terms of solvation properties, in the absence of po- 
larizability change of the solute, AAo/2 should be equal 
to A = Ai = A2. The first lines in Tables HI and HT1 corre- 
spond to exactly this situation. The difference between 
two reorganization energies, as well as deviations of both 
of them from AAo/2 is another indication of the non- 
Gaussian statistics of the solvent polarization leading to 
non-linear solvation. Despite these complications, the di- 
rect calculations of the Stokes shift from simulations com- 
pare semi-quantitatively to the results of applying ea 1381 
(Tablcs[l]and[nJ). The analytical model developed in Sec- 
tion therefore, captures the basic thermodynamics 
of ET with markedly different reorganization energies. 

TABLE II: Energetic parameters (eV) of the charge recombi- 
nation reaction (Figure |Sj for which 002 is varied and a g — 
is kept constant; Pqo/a = 10 in the charge-separated state. 

<W<t 3 \ 3 a A 4 AX " AAV 



0.0 
0.2 
0.3 
0.4 
0.5 
0.6 
0.7 
0.8 



A3 

0.369 0.279 

0.654 0.284 

1.115 0.287 

1.466 0.290 

2.270 0.293 

2.562 0.296 

3.240 0.299 

7.466 0.303 



0.649 
0.911 
1.115 
1.466 
2.270 
2.562 
3.240 
4.079 



0.642 
0.871 
1.051 
1.352 
1.728 
2.357 
2.928 
3.446 



"Calculated from the simulation data assuming /3 = 40 eV 
''From th e sim ulation data. 
c From eg 1551 



IV. DISCUSSION 

The basic design of an artificial photosynthetic device, 
as advances by Meyer and co-workers^2^i is shown in 
Figure El It anticipates the creation of ordered ar- 
rays where donors and acceptors are connected to cat- 
alytic sites where high-energy reactions can occur (e.g., 
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splitting of water or reduction of carbon dioxide to car- 
bohydrates). This design requires efficient separation 
of the electron and the hole. Unless this separation 
is achieved by high mobility of each charge in molec- 
ular arrays^ or in valence and conduction bands of a 
semiconductor^ high branching ratio between charge 
separation and charge recombination is required at each 
step of productive charge transfer to facilitate the redox 
reactions which are normally significantly slower than ET 
steps. 



Red, 




FIG. 11: Reaction center model^^i The unit C is photoex- 
cited to obtain the electron from donor D and transfer it to 
acceptor A. The recombination rates k r at each stage should 
be sufficiently low to allow normally slow catalytic reduction 
(Oxi — >Redi) and oxidation (Red2 -^0x2) reactions to occur. 

The model presented here offers some novel opportu- 
nities in terms of varying the parameters of ET reactions 
and activation barriers. The model emphasizes the cou- 
pling of the solute polarizability change, achieved in the 
course of electronic transition, to the electric field at the 
position of the donor-acceptor complex. The electric field 
has two components, the reaction field of the polar envi- 
ronment Ri and the macroscopic field F m . The coupling 
of the overall field Ri + F m to the polarizability change 
creates the effective dipole moment change of the donor- 
acceptor complex fea 1140 




FIG. 12: Free energy surfaces of photoinduced ET in the 
Marcus-Hush picture (a) and in the present model with dif- 
ferent reorganization energies (b). The free energy surfaces 
are plotted against the reaction coordinate corresponding to 
the energy gap between the charge-separated and photoex- 
cited states. The parameters are chosen to show the activa- 
tionless pathway from the photoexcited state 1 to the charge- 
separated state 2 and from 2 to the ground state 3. The 
free energy surfaces 2 and 3 coincide in the Marcus-Hush pic- 
ture, but are distinct in the present model. In (a) the dashed 
line indicates the free energy surface for the normal-region 
charge recombination when energetic efficiency is below 50%; 
hv indicates the energy of photoexcitation. The reorganiza- 
tion energies are: A = 1 eV (a) and Ai = 0.2 eV, A2 = 2.0 eV, 
A 3 = 1.0 eV, and A 4 = 0.5 eV (b). 



Am e ff,i = Am + Aa(Ri + F m ) (41) 

The reaction field Ri, which depends on the electronic 
state of the donor-acceptor complex, is responsible for 
the distinction between the forward and backward reor- 
ganization energies, Ai ^ A2 and A3 7^ A4 as is shown in 
Figure E3 and listed in Tables [Q and HJ Note that the 
asymmetry in the reorganization energies up to a factor 
of 25 obtained in our simulations ( Table HT|) is the largest 
ever observed in either computer or laboratory experi- 
ment. 

A note on the dual nature of the reorganization energy 
is relevant here. The reorganization energy for each ET 
state is calculated on the configurations of the solvent 
in equilibrium with the donor-acceptor complex in that 
given state. Therefore, the reorganization energy car- 
ries information about a particular electronic state of the 
donor-acceptor complex. On the other hand, the value 
that is being averaged is the change in the interaction 
potential, which is a characteristic of a given transition. 
As a result of this duality, the electronic charge-separated 



state of the donor-acceptor complex D + -A~ can be char- 
acterized by two, generally unequal, reorganization en- 
ergies A2 and A3 reflecting transitions to two different 
electronic states, photoexcited and ground. 

In the present picture, all reorganization energies Ai, 
A2, A3, and A4 are allowed to be different, producing 
reach energetics of photoinduced charge separation and 
charge recombination. The difference in scenarios for 
photoinduced ET in the Marcus-Hush model and in the 
present formulation is illustrated in Figure 1121 The 
plots show the free energy surfaces vs the reaction co- 
ordinate corresponding to the energy gap X between 
the charge-separated and photoexcited states. The en- 
ergy gap for the charge recombination reaction is then 
—AFqs — AFcr + X, where AFcs.cr are the free energy 
gaps for charge-separation and charge-recombination re- 
actions. In the Marcus-Hush model, activationless charge 
separation and charge recombination is achieved at the 
photoexcitation energy hv = 2 A (Figure IT2T a)h There- 
fore, the desire to move the recombination reaction to the 
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normal region requires energetic efficiency of the photo- 
synthetic apparatus to be below 50% (dashed line in Fig- 
ure II 2f a)). This situation changes when different reorga- 
nization energies are allowed in the construction of the 
free energy surfaces. The plots in Figure IT^T b) are made 
with A2/A1 = 10 and A3/A4 = 5. When activationless re- 
actions are realized for 1 — > 2 and 2 — > 3 transitions, the 
asymmetry of the energy gap law leads to a higher ener- 
getic efficiency of 70%. The disadvantage of this scheme 
is a shallow shape of the charge-separated free energy sur- 
face making activation energy for charge recombination 
rather weakly dependent on changes in the free energy 
gap. 

The role of the reaction field factor in eq^Jwill dimin- 
ish in weakly polar media characteristic of protein elec- 
tron transfer. The effect of a strong local electric field 
F m then becomes more important. If the term AaRi 
is small compared to Am, we can simplify the effective 
dipole moment change to the form 

Am cff = Am + AaF m (42) 

Since Am c ff does not depend on the electronic state, the 
forward and backward reorganization energies are the 
same (Ai = A2, A3 = A4), according to the Marcus-Hush 
theory. However, the factor AaF m can be responsible 
for the difference in the reorganization energies between 
charge separation and charge recombination reactions. 




-3-2-10123 
X/eV 



FIG. 13: Free energy surfaces in the Marcus-Hush picture 
when the reorganization energy for charge recombination 
A3 = A4 = 2.0 eV is higher than reorganization energy for 
charge separation Ai = A2 = 0.5 eV due to the coupling of 
the polarizability change to the macroscopic electric field (eqs 
d and |13 . 



Imagine a situation in which the photoexcited state 
is significantly more polarizable than the ground state 
and the charge-separated state has about the same po- 
larizability as the photoexcited state, Aacs — 0- Then, 
in weakly polar solvents, Ai ~ A2 and both reorganiza- 
tion energies are small because of the low polarity of the 
medium. For the charge recombination reaction, Am < 
and Aq;cr < 0. The permanent dipole Am and the 
induced dipole AaF m in eq 021 add up when F m > 
and subtract when F m < 0. In the former case, the 
charge-recombination reorganization energy A3 ~ A4 is 
higher than the charge-separation reorganization energy 
Ai ~ A2. Note that this picture does not contradict to the 
energy conservation requirement given by eq^ since the 
photoexcited state D*-A and the ground state D-A are 
two different electronic states characterized by different 
polarizabilitics. 

The normal-region ET can be realized within the 
Marcus-Hush picture when the reorganization energy is 
greater than the charge-recombination free energy gap 
A 3 > AFcr (Figure IT3|) . As in Figure Eta), the ener- 
getic efficiency is about 50%, and it drops when charge 
recombination is shifted into the normal region. How- 
ever, even in this limit of reduced flexibility of altering 
the ET parameters, the presence of the macroscopic elec- 
tric field in the expression for the reorganization energy 
fens HI and H2|) opens the door to manipulations of reor- 
ganization parameters through proper design of mutual 
signs of the polarizability change and the direction of the 
macroscopic field. 

We can conclude that the combination of a highly po- 
larizable donor-acceptor complexe with macroscopic elec- 
tric field creates a principal possibility for suppressing 
the recombination reaction in photosynthctic ET. From a 
more general perspective, the notion of non-equal reorga- 
nization energies and non-parabolic free energy surfaces, 
compliant with the condition of energy conservation (eq 
, opens the door to new models of ET activation which 
may eventually lead to a practical solution of the photo- 
synthesis problem. 
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